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I define a set of wavefunctions for SU(A'') lattice antiferromagnets, analogous to the valence bond 
solid states of Affieck et al. [J] , in which the singlets are extended over A^-site simplices. As with the 
valence bond solids, the new simplex solid (SS) states are extinguished by certain local projection 
operators, allowing one to construct Hamiltonians with local interactions which render the SS states 
exact ground states. Using a coherent state representation, I show that the quantum correlations in 
each SS state are calculable as the finite temperature correlations of an associated classical model, 
with A-spin interactions, on the same lattice. In three and higher dimensions, the SS states can 
spontaneously break SU(A) and exhibit A-sublattice long-ranged order, as a function of a discrete 
parameter which fixes the local representation of SU(A). I analyze this transition using a classical 
mean field approach. For N > 2 the ordered state is selected via an 'order by disorder' mechanism. 
As in the AKLT case, the bulk representations fractionalize at an edge, and the ground state entropy 
is proportional to the volume of the boundary. 

PACS numbers: 75.10.Hk, 75.10.Jm 



I. INTRODUCTION 

At the classical level, the thermodynamic properties 
of fcrromagnets and antiferromagnets are quite similar. 
Both states break certain internal symmetries, whether 
they be discrete or continuous, and often crystalline point 
group symmetries as well. Antiferromagnetism holds the 
interesting possibility of frustration, which can lead to 
complex behavior even at the classical level. 

Quantum mechanics further distinguishes antiferro- 
magnetism as the more interesting of the two phenom- 
ena. Quantum fluctuations compete with classical order- 
ing, and many models of quantum antiferromagnetism 
remain disordered even in their ground states. The rea- 
son is that on the local level, quantum antiferromagnets 
prefer distinctly non-classical correlations, in that they 
form singlets, which are superpositions of classical states. 
For a S" = i Heisenberg antiferromagnet on a bipartite 
lattice, theorems by Marshall and by Lieb and Mattis 
[1] rigorously prove that the ground state is a total spin 
singlet: S^^^ = 0. Any total singlet can be expanded in a 
(nonorthogonal) basis of valence bonds, which are singlet 

pairs {ij) = 2^"^/^(| Tj ij ) ^ Uitj)) extending between 
sites i and j. The most probable singlets are between 
nearest neighbors, which thereby take full advantage of 
the Heisenberg interaction J ■ Sj and achieve a mini- 
mum possible energy Eq = — | J for that particular link. 
Taking linear combinations of such states lowers the en- 
ergy, via dclocalization, with respect to any fixed singlet 
configuration; this is the basic idea behind Anderson's 
celebrated resonating valence bond (RVB) picture If 
one allows the singlet bonds to be long-ranged, such a 
state can even possess classical Neel order 

Taking advantage of quantum singlets, one can 
construct correlated quantum-disordered wavefunctions 
which are eigenstates of local projection operators. This 
feature allows one to construct a many-body Hamiltonian 



which renders the parent wavefunction an exact ground 
state, typically with a gap to low-energy excitations. Per- 
haps the simplest example is the Majumdar- Ghosh (MG) 
model for a. S — spin-chain Q, the parent state of 
which is given by alternating singlet bonds, viz. 



(1) 



The key feature to | 5*) is that any consecutive trio of 
sites (n, n -I- 1, ?7 + 2) can only be in a state of total spin 
S ~ ^ ~ there is no 5 = | component. Thus, | ^'^ is an 
eigenstate of the projection operator 



3/2 



(n, n-|-l,n + 2) 



(2) 



with zero eigenvalue, and an exact ground state for 
n = JY.n^3/2('^^'^ + 1'" + 2)- As l^f) breaks lat- 
tice translation symmetry, a second (degenerate) ground 
state follows by shifting | by one lattice spacing. Ex- 
tensions of the MG model to higher dimensions and to 
higher spin, where the ground state is again of the Kekule 
form, i.e. a product of local valence bond singlets, were 
discussed by Klein Q. 

Another example is furnished by the valence bond 
solid (VBS) states of Affleck, Kennedy, Lieb, and Tasaki 
(AKLT) The general AKLT state is compactly writ- 
ten in terms of Schwinger boson operators 
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which assigns m singlet creation operators to each link 
of a lattice C. The total boson occupancy on each site is 
zM , where z is the lattice coordination number; in the 
Schwinger representation this corresponds to 25". Thus, a 
discrete family of AKLT states with S = \zM is defined 
on each lattice, where M is any integer. The maximum 
total spin on any link is then S"!^^^ = 2S — M , and any 
Hamiltonian constructed out of link projectors for total 



2 



spin S G [25 ~ M + 1 , 28] , with positive coefRcients, 
renders | Af)) an exact zero energy ground state. 

The elementary excitations in these states were treated 
using a single mode approximation (SMA) in ref. 0. 

The ability of two spins to form a singlet state is a spe- 
cial property of the group SU(2). Decomposing the prod- 
uct of two spin-S* representations yields the well-known 
result, 



S (S) S = 0®1®2< 



25, 



(4) 



vacua. More recently, Greiter and Rachel considered 
SU(iV) valence bond spin chains in both the fundamental 
and other representations, constructing their correspond- 
ing Hamiltonians, and discussing soliton excitations. Ex- 
tensions of Klein models, with Kckule ground states con- 
sisting of products of local SU(A^) singlets, were discussed 
by Shen |13| . and more recently by Nussinov and Ortiz 
[l4| . An SU(4) model on a two leg ladder with with a 
doubly degenerate Majumdar- Ghosh type ground state 
has been discussed by Chen et al. [l^ . 



and there is always a singlet available. If we replace 
SU(2) by SU(3), this is no longer the case. The repre- 
sentations of SU(N) are classified by {N — l)-row Young 
tableaux {li,l2, ■■■ with Ij boxes in row j, and 

with ^1 > ^2 ^ ■ ■ ■ ^ ^N-i- "^^^ product of two funda- 
mental (1,0) representations of SU(3) is 



□ ® □ 
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(5) 



which docs not contain a singlet. One way to rescue the 
two-site singlet, for general SU(7V), is to take the product 
of the fundamental representation N with the antifun- 
damental N. This yields a singlet plus the (iV^ — 1)- 
dimcnsional adjoint representation. In this manner, gen- 
eralizations of the SU(2) antiferromagnet can be defined 
in such a manner that the two-site valence bond structure 
is preserved, but only on bipartite lattices [1, [l^l- 

Another approach is to keep the same representation 
of SU(iV) on each site, but to create SU(A^) singlets ex- 
tending over multiple sites. When each site is in the 
fundamental representation, one creates A'^-site singlets. 



(%) 1 0) , 



(6) 



where 6jj(i) creates a Schwinger boson of flavor index a 
on site i. The SU(A'^) spin operators may be written in 
terms of the Schwinger bosons as 



(7) 



with Tr (5) = 0, for the general symmetric (p, 0) repre- 
sentation. These satisfy the SU(A^) commutation rela- 
tions. 



(8) 



Extended valence bond solid (XVBS) states were first 
discussed by Affleck et al. in ref. [11|. In that work, 
SU(2A^) states where N — mz were defined on lattices of 
coordination number z, with singlets extending over z + l 
sites. Like the MG model, the XVBS states break lat- 
tice translation symmetry t and their ground states are 
doubly degenerate; they also break a charge conjugation 
symmetry C, preserving the product tC. In addition to 
SMA magnons, the XVBS states were found to exhibit 
soliton excitations interpolating between the degenerate 



Shen also discussed a generalization of Anderson's 
RVB state to SUfiV) spins, as a prototype of a spin- 
orbit liquid state ^M,- A more clearly defined and well- 
analyzed model was recently put forward by Pankov, 
Moessner, and Sondhi [l^j who generalized the Rokhsar- 
Kivelson quantum dimer model p7| to a model of res- 
onating singlet valence plaquettes. Their plaquettes are 
Af-site SU(A^) singlets (A^ = 3 and A^ = 4 models 
were considered), which resonate under the action of the 
SU(A'^) antiferromagnctic Heisenberg Hamiltonian, pro- 
jected to the valence plaquette subspacc. The models and 
states considered here do not exhibit this phenomenon 
of resonance. Rather, they are described by static "sin- 
glet valence simplex" configurations. Consequently, their 
physics is quite different, and in fact simpler. For exam- 
ple, with resonating valence bonds or plaquettes, one can 
introduce vison excitations [isj which are Z2 vortex ex- 
citations, changing the sign of the bonds or plaquettes 
which are crossed by the vortex string (l9| . For simplex 
(or plaquette) solids, there is no resonance, and the vison 
does not create a distinct quantum state. The absence of 
'topological quantum order' in Klein and AKLT models 
has been addressed by Nussinov and Ortiz . 

Here I shall explore further generalizations of the 
AKLT scheme, describing a family of 'simplex solid' (SS) 
states on A^-partite lattices. While the general AKLT 
state is written as a product over the links of a lattice £, 
with M singlet creation operators applied to a given link, 
the SS states, mutatis mutandis^ apply M SU(A'^) singlet 
operators on each A^-simplex. Each site then contains 
an SU(A'^) spin whose representation is determined by 
M and the lattice coordination. Furthermore, as is the 
case with the AKLT states, the SS states admit a simple 
coherent state description in terms of classical CP^~^ 
vectors. Their equal-time quantum correlations are then 
computable as the finite temperature correlations of an 
associated classical model on the same lattice. A classical 
ordering transition in this model corresponds to a zero- 
temperature quantum critical point as a function of M 
(which is, however, a discrete parameter). I argue that 
the ordered SS states select a particular ordered struc- 
ture via an 'order by disorder' mechanism. Finally, I 
discuss what happens to these states at an edge, where 
the bulk SU(A^) representation is effectively fractional- 
ized, and a residual entropy proportional to the volume 
of the boundary arises. 



II. SIMPLEX SOLIDS 

Consider an iV-site simplex 7^, and define the SU(-/V) 
singlet creation operator 



n 



"^1 (A)' 



(9) 



where i = 1, . . . , iV labels the sites on the simplex. 
Any permutation ^ of the labels has the trivial conse- 
quence of TZ^p sgn(/^) 7?.|n. Next, partition a lattice C 
into A^-site simplices, i.e. into N sublattices, and define 
the state 



*(£;M))=^(7^^^|0) 



(10) 





FIG. 1: Top: SU(3) SS state on two-leg zigzag chain. Each 
site is in the 10-dimensional totally symmetric representation 
with 3 boxes (the red dots). Bottom: SU(4) SS state on two- 
leg ladder of tetrahedra. Each site is in the 10-dimensional 
totally symmetric representation with 2 boxes. 



where M is an integer. Since each 7i}p operator adds one 
Schwinger boson to every site in the simplex, the total 
boson occupancy of any given site is p = CM, where 
C is the number of simplices associated with each site. 
For lattices such as the Kagome and pyrochlore systems, 
where two neighboring simplices share a single site, we 
have C = 2. For the tripartite triangular lattice, C = 
3. Recall that each site is in the (p, 0) representation of 
SU(A^), with one row of p boxes. 

Two one-dimcnsional examples are depicted in fig. [TJ 
The first is defined on a two-leg zigzag chain. The chain is 
partitioned into triangles as shown, with each site being 
a member of three triangles. Each triangle represents 
a simplex F and accommodates one power of the SU(3) 
singlet creation operator TZ]^. Thus, for this state we have 
iV = 3 and C = 3. With M = 1 then, the local SU(3) 
representation on each site is (3,0), i.e. mm . which is 
10-dimcnsional. The M = 1 case is in fact a redrawn 
version of the state defined by Greiter and Rachel in cqn. 
52 of ref. [l2|- For this state, a given link may be in any 
of four representations of SU(3), or a hnear combination 
thererof: 



10 



(11) 
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35 



28 



Thus, the zigzag chain SU(3) SS state is an exact zero 
energy eigenstate of any Hamiltonian of the form 







(12) 



p (^.^) + J3 P| 1 1 1 1 i (^i) 



with positive coefficients J^, Jj, and Jg. 

The SU(4) SS chain in fig. [T]is topologically equivalent 
to a chain of tetrahedra, each joined to the next along 
an opposite side. Thus, C = 2 and for the M = 1 parent 



state, each site is in the 10-dimensional n~l representa- 
tion. From 



m 
10 



20 



F 

45 



we can construct a Hamiltonian, 



(13) 



(14) 



E "^2 Prrrn(u) + ^3 P^(*i) > 



rungs 



again with positive coefficients J^, J2, Jg and J4, which 
renders the wavefunction | ^E*^ an exact zero energy 
ground state. For both this and the previously discussed 
SU(3) chain, the ground state is nondegenerate. 



A. SU(iV) Casimirs 

For a collection of K spins, each in the fundamental of 
SU(7V), we write 



K 



E^^(^) 



(15) 



fc=i 



From the spin operators S^, one can construct A — 1 
Casimirs, C^"' = ^Tr(5"), with n = 2,...,A. The 
eigenvalues of C^"-* for totally symmetric (+) and totally 
antisymmetric (— ) representations of p boxes were ob- 
tained by Kobayashi [20|: 



((A^Tl)(^±p) 



(16) 



The Casimirs can be used to construct the projectors 
onto a given representation as a polynomial function of 



4 



the spin operators. In order to do so, though, we will need 
the eigenvalues for all the representations which occur in 
a given product. Consider, for example, the case of three 
SU(3) objects, each in their fundamental representation. 
We then have 



□ ® □ ® □ 



(17) 



The eigenvalues of the quadratic and cubic Casimirs are 
found to be 



Therefore 







2 M2) 

3 ^ 
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P^i^jk) = -2 + 



ran 



fe-fc) = l-iC(2) + iC(3) 



(18) 
(19) 
(20) 



Expressing the projector P ^-^ j^ijk) in terms of the local 
spin operators 5*^(0; I find 

P^(u/c) = 1 - i Tr (4,) + Tr (4,) (21) 
= - i Tr [S{i) S{j) + S{j) S{k) + S{k) S{i)] 
+ i Tr [Sit) Sij) S{k) + S{k) S{j) S{z)] - ^ . 

The projector thus contains both bilinear and trilinear 
terms in the local spin operators. One could also write 
the projector in terms of the quadratic Casimir C*-^-* only, 
as 



(22) 



This, however, would result in interaction terms such 
as Ti-[S{i)S{j)] ■ TT[S{j)S{k)], which is apparently 
quadratic in S{j). For the □ representation, however, 
products such as Sp{i)Sf^{i) can be reduced to hnear 
combinations of the spin operators S!^{i), as is familiar 
in the case of SU(2). This simplification would then re- 
cover the expression in eqn. 1211 



III. SS STATES IN d > 2 DIMENSIONS 

A. Kagome Lattice 

The Kagome lattice, depicted in fig. [21 is a two- 
dimensional network of corner-sharing triangles, with 
( = 2. It naturally accommodates a set of -/V = 3 SS 
states. The simplest example consists of SU(3) objects 
in the fundamental representation at each site, and places 




FIG. 2: SU(3) simplex solid states on the Kagome lattice. 
Applying the singlet operator TZ\, to all the up (down) tri- 



angles generates the state N/ 



Al V 



Applying TZ]^ to all the 



triangular simplices generates the state | of eqn. [2S] The 
twofold coordinated yellow sites at the top form a (10) edge 
(see sec. IVIII|I . 



SU(3) singlets on all the upward-pointing triangles (see 
fig.©: 



*A)=n^Aio> • 



(23) 



A 



The lattice inversion operator T then generates a degen- 
erate mate, | = T- \ 4*^) Both states are exact zero 



energy eigenstates of the Hamiltonian 



n 



= E 

(y"fe)6120° 



(24) 



where the sum is over all 120° trios (ijk); there are six 
such (ijk) trios for every hexagon. Since two of the three 
sites in each trio are antisymmetrized. The fully symmet- 
ric I I I I representation is completely absent. This model 
bears obvious similarities to the MG model: its ground 
state is a product over independent local singlets, hence 
there are no correlations beyond a single simplex, and it 
spontaneously breaks a discrete lattice symmetry (in this 
caseX) [U. 

If we let the singlet creation operators act on both the 
up- and down-pointing triangles, we obtain a state which 
breaks no discrete lattice symmetries. 



^)=n^AiKio) 



(25) 
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For this state, each site is in the 6-dimensional i i i rep- 
resentation. On any given link, then, there are the fol- 
lowing possibilities: 



(26) 



The fact that each link belongs to either an A or v 
simplex, and the fact that a singlet operator TZ 



t 

A/V 



IS 



associated with each simplex, means that no hnk can be 
in the fuUy symmetric i i i i i representation. Thus, | 4*) is 
an exact, zero-energy eigenstate of the Hamiltonian 



(27) 



The states | ^'^), | ^',^), and | arc depicted in fig. [3 
The actions of the quadratic and cubic Casimirs on the 
possible representations for a given link arc given in the 
following table: 



Rep" 


C7(2) 


C(3) 




10/3 


-80/27 




16/3 


64/27 


1 1 1 1 1 


28/3 


352/27 



Note that 



C(3) 



(28) 



hence the two Casimirs are not independent here. We 
can, however, write the desired projector, 



rrm 



(u) 



^(2 



24 



Tr[5(j)5(j)]) +^TT[S{i)Sij)] 



(29) 



_5_ 

27 ' 



as a bilinear plus biquadratic interaction between neigh- 
boring spins. To derive this result, we write C*-^^ = 
^Ti- [S{i) + S{j)]^, whence 

= Tr [S{i) Si,])] + Tr (S^) 

^Tr[S{t)Sij)] +p{N + p- 



' N 



(30) 



Next, consider the pyrochlore lattice in fig. [31 This 
lattice consists of corner-sharing tetrahedra, with C = 2, 
and naturally accommodates an = 4 SS state of the 
form 



*) = n 0) 

tetrahedra 



(31) 



Like the uniform SU(3) SS state on the Kagomc lattice, 
this SU(4) state describes a lattice of spins which are in 
the representation on each site; in the SU(4) case 

this representation is 10-dimensional. From eqn. 1131 we 
see that each link, the sites of which appear in some sim- 
plex singlet creation operator, cannot have any weight 




FIG. 3: The pyrochlore lattice and its quadripartite structure. 

in the 35-dimcnsional totally symmetric Mill repre- 
sentation. Hence, once again, the desired Hamiltonian is 
that of eqn. [271 For SU(4), 

C(2)(g)=6 , C(2)(^)=8 , C(2)(™) = 12, 
and so 



rrm 



(^,) = ^(C(^)-6)(C(2)-8) 



(32) 

= ^ (Tr [Siz) S{j)\)\ I Tr [S{z) S{j)\ + \ , 

Indeed, there is a rather direct correspondence between 
the possible SU(3) SS states on the Kagome lattice and 
the SU(4) SS states on the pyrochlore lattice. For exam- 
ple, one can construct a model with a doubly degenerate 
ground state, similar to the MG model, by associating 
the simplex singlet operators Ti)p with only the tetrahe- 
dra which point along the (111) lattice direction. 

Finally, consider SU(4) states on the square lattice, 
again in the □□ representation on each site. We can once 
again identify the exact ground state of the ^^jjjjij) 
Hamiltonian of eqn. [23 In this case, the ground state 
is doubly degenerate, and is described by the 'planar py- 
rochlore' configuration shown in fig. [31 



IV. MAPPING TO A CLASSICAL MODEL 

The correlations in the SS states are calculable using 
the coherent state representation. From results derived 
in the Appendix I, the coherent state SS wavefunction is 
given by 



^[{-zm]=cX{[RMr,),...rz{rj,)) 

r 

where C is a normalization constant, and 



M 



(33) 



(34) 
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As shown in Appendix I, the matrix clement 



I '^K I V') operator 



■b 



t'jv 

AT 



(39) 




FIG. 4: One of two doubly degenerate ground states for 
the Hamiltonian of eqn. [27] apphed to the square lattice, 
where each site is in the (2,0,0) representation of SU(4). 
The squares with crosses (or, equivalently, tetrahedra) repre- 
sent singlet operators e"'^''^ ba{i) b^ij{j) btf{k) b\{l) on the sim- 
plex (ijkl). The resulting 'planar pyrochlore' configuration is 
equivalent to a checkerboard lattice of tetrahedra. 



Here I have labeled the N sites on each simplex F by an 
index i running from 1 to N. 

Note that the coherent state probability density is 



\^\' = \C\'l[\R, 



|2M 



(35) 



may be computed as an integral with respect to the mea- 
sure dfi (on each site) of the product (j>{z) 'ip{z) of coherent 
state wavefunctions multiplied by the kernel 



TK{{K}Ahi}) 



(N -1 



■P- 



K)\ 



(40) 

Thus, the quantum mechanical expectation values of Her- 
mitian observables in the SS states are expressible as 
thermal averages over the corresponding classical Hamil- 
tonian H^^ of eqn. [37l The SS and VBS states thus 
share the special property that their equal time quan- 
tum correlations are equivalent to thermal correlations 
of an associated classical model on the same lattice, i.e. 
in the same number of dimensions. 

In this paper I will be content to merely elucidate 
the correspondence between quantum correlations in 
I 4'(£ ; M)) and classical correlations in H^^. An applica- 
tion of this correspondence to a Monte Carlo evaluation 
of the classical correlations will be deferred to a future 
publication. 



and that 



".•-„,ft-/3„Q^^^^(rj...Q^ , {r^) (36) 



where Q^pii) = z^(i)zp{i). Writing \^\'^ = e^^^^^'^, we 
see that that the probability density may be written as 
the classical Boltzmann weight for a system described by 
the classical Hamiltonian 



= -^ln|i?. 



(37) 



V. SINGLE MODE APPROXIMATION FOR 
ADJOINT EXCITATIONS 

Following the treatment in ref. [1], I construct trial 
excited states at wavevector k in the following manner. 
First, define the operator 

(41) 

i 

R 



at a temperature T = 1/M [1, [24I. The classical inter- 
actions are iV-body interaction, involving the matrices 
Q^p{i) on all the sites of a given A^-site simplex, summed 
over all distinct simplices. For N = 2, this results in a 
classical nearest- neighbor quantum antiferromagnet 3, 
with 



where i? is a Bravais lattice site and i labels the basis 
elements. Here, rj^ is for the moment an arbitrary set of 
complex- valued parameters and is the total number of 
unit cells in the lattice. The operators (f'api^) transform 
according to the {N'^ — l)-dimensional adjoint represen- 
tation of SU(A^). Next, construct the trial state. 



1 



(38) 



with n = z'crz, where cr is the vector of Pauli matrices. 
This general feature of pair product wavefunctions of the 
Bijl-Feynman, Laughlin, and AKLT form is thus valid for 
the SS states as well. 



|</>)^,/.„,(fe)|*) 



(43) 



and evaluate the expectation value of the Hamiltonian in 
this state: 



^SMA(fc) 





H 









(44) 
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Here, f^jik) and s^j{k) are, respectively, the oscillator 
strength and structure factor, given by 

= K*! K./jC^) ' [H, I*) (45) 

s,^.(fc) = (vl/|4,^(fc)0^.„^(fc)|vl/) . (46) 

Here I have assumed that i7 is a sum of local projec- 
tors, and that -ffj^') = 0. Treating the r]^ parameters 
variationally, one obtains the equation 



4(fc)'?,— ^SMA(fc)s,,(fc)^?, 



(47) 



The lowest eigenvalue of this equation provides a rigorous 
upper bound to the lowest excitation energy at wavevec- 
tor k. The result is exact if all the oscillator strength 
is saturated by a single mode, whence the SMA label. 
When is quantum-disordered, the SMA spectrum is 
gapped. When j^*) develops long-ranged order (for suf- 
ficiently large M parameter, and in d > 2 dimensions), 
the SMA spectrum is gapless. 



VI. MEAN FIELD TREATMENT OF 
QUANTUM PHASE TRANSITION 

The classical Hamiltonian of eqn [37] exhibits a global 
SU(iV) symmetry, where z^{i) — > U^^, z^,(i) for every 
lattice site i. Since the interactions are short-ranged, 
there can be no spontaneous breaking of this symmetry 
in dimensions d < 2. In higher dimensions, the classical 
model can order at finite temperature, corresponding to 
a quantum ordering at a finite value of m. For the AKLT 
states, where N = 2, this phase transition was first dis- 
cussed in ref . Q . I first discuss the N ~2 case and then 
generalize to arbitrary N > 2. 



A. iV = 2 : VBS States 

Consider the N = 2 case, which on a lattice of coor- 
dination number z yields a family of wavefunctions de- 
scribing S = \Mz objects with antifcrromagnctic corre- 
lations. These are the AKLT valence bond solid (VBS) 
states. We have 



3in(§) e*"^ 




(48) 

where n = z^ crz is a real unit vector, (cr arc the Pauli 
matrices). Since 

e".ft Q (^j Q (^j = 1(1 - . O , (49) 



the effective Hamiltonian is 



El- 



1 — n,, • h„ 



(50) 



The sum is over all links on the lattice. I assume the 
lattice is bipartite, so each link connects sites on the A 
and B sublattices. I now make the mean field Ansatz 

fij^ = m + Shj^ , fi-g = m + Sfi^ (51) 

and expand H in powers of 5h^. Expanding to lowest 
nontrivial order in the fluctuations Sn^, we obtain a mean 
field Hamiltonian 

H"- ^E,-KY.f<'^^^-T.^3 ^ (52) 



where is a constant and 



(53) 



is the mean field. Here z is the lattice coordination num- 
ber. The self-consistency equation is then 



m = 



= jdhh e^^ / J dn e^A '^/^ , (54) 



which yields 



ctnh A 



1 
A 



T 1 



(55) 



The classical transition occurs at T^"^^ = ^z, so the SS 
state exhibits a quantum phase transition at il/"'' ~ 
3z~^. For M > the SS exhibits long-ranged two- 
sublattice Neel order. On the cubic lattice, the mean 
field value A/"'' = i suggests that all the square lattice 
VBS states, for which the minimal spin is S* = 3 (with 
M = 1), are Neel ordered. Since the mean field treat- 
ment overestimates due to its neglect of fluctuations, 
I conclude that the true is somewhat greater than 
3z^^, which leaves open the possibility that the minimal 
M = 1 VBS state on the cubic lattice is a quantum dis- 
ordered state. Whether this is in fact the case could be 
addressed by a classical Monte Carlo simulation. 



B. iV > 2 : SS States 

For general N, I write 



(56) 



where the average is taken with respect to = 
e~^t:i/^. To maximize jvE'P, i.e. to minimize H^^, choose 
a set {PJ^^} of N mutually orthogonal projectors, with 
CT = 1, . . . , A^. The projectors satisfy the relations 



p-^ P'^ = S, P-^ , 
and can each be written as 



(57) 
(58) 
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where is a set of TV mutually orthogonal CP^~^ 

vectors. Then if z{r^) ~ oj'^ for each site in the 
simplex, we have Rp = e"* , where rj is an arbitrary phase, 
and \Rp\'^ = 1. One then writes 

2a/3^(Qa/3(^J> = ^^./3 + "^(Pa/3-^'5„/3) • (59) 

Here m £ [0, 1] is the order parameter, analogous to the 
magnetization. When m = 0, no special subspacc is se- 
lected, and the correlations are isotropic. When m = 1, 
the Q-matrix is a projector onto the one-dimensional sub- 
space defined by uj'^. Note that {Tr Q{r^)) ~ 1. as it 
must be. 

At this point, there remains a freedom in assigning the 
vectors {io'^} to the sites {F^} of each simplex. Consider, 
for example, the iV = 3 case on the Kagome or triangu- 
lar lattice. The lattice is tripartite, so every A sublattice 
site has 2 (Kagome) or 3 (triangular) nearest neighbors 
on each of the B and C sublattices. However, as is well- 
known, the individual sublattices may have lower transla- 
tional symmetry than the underlying triangular Bravais 
lattice. Indeed, the sublattices may be translationally 
disordered. I shall return to this point below. For the 
moment it is convenient to think in terms of N sublat- 
tices each of which has the same discrete symmetries as 
the underlying Bravais lattice. 

Expanding H^^ to lowest order in the fluctuations 
5Q^Ji) on each site, and dropping terms of order 



and higher, I obtain the mean field Hamiltonian, 



(60) 



where Eq is a constant. The mean field h^^ii) is site- 
dependent. On a site, we have 

(61) 

- i^) Lfi + (m) P^„) / Rj, (m) , (62) 
where 

^^(m) = (A^-2)! ^i^--^,„WL^ 



In Appendix II I show that 



3=0 



N-2 



B^{m)^{N-2)\ J2 



J! 



2^ + 1 /I 



3=0 



j\ \ N 



N 



N-3-2 



^ 1/1 \N~3 



3=0 



(63) 
(64) 

(65) 



Note that ^^^(0) = 0, B^(l) = 1, and R{1) = 1. 
The mean field Hamiltonian is then 



i 



(66) 



where h^p[i) = h'^^^\ where a{i) labels the projector 
associated with site i. The self-consistency relation is 
obtained by evaluating the thermal average of Qapii)- 
With = Iz^p, I obtain 



1 A^- 1 



N 



N 



J dx X (1 - x)^-2 exp(C 6^(m) x/T) 



1 ' 
jdx (1 - x)^-2 exp(C 6^(m) x/T) 



(67) 



where hj^{m) = Bj^{m)/Rj^{m). It is simple to see that 
TO = is a solution to this mean field equation. To find 
the critical temperature T^, expand the right hand side 
in powers of to for small m. To lowest order, one finds 



N 



N -I 



(68) 



The value of T^^ is determined by equating the coeffi- 
cients of TO on either side of the equation. I find 



M^^iN, C) 



1 



c 



(69) 



This agrees with our previous result M"^ = i for the 
N ~ 2 state on the cubic lattice, for which C 6. The 
A^ = 3 SS on the Kagome lattice cannot develop long- 
ranged order which spontaneously breaks SU(A'^), owing 
to the Mermin- Wagner theorem. For the A^ = 4 SS on 
the pyrochlore lattice, our mean field theory analysis sug- 
gests that the SS states are quantum disordered up to 
A/ sa ^ . Note that expression for M"'' reflects a compe- 
tition between fluctuation effects, which favor disorder, 
and the coordination number, which favors order. The 
numerator, A^^ — 1, is essentially the number of directions 
in which Q can fluctuate about its average Q; this is the 
dimension of the Lie algebra su(A^). 



VII. ORDER BY DISORDER 



At zero temperature, the classical model of eqn. [37] is 
solved by maximizing \Rp\ for each simplex F. This is 
accomplished by partitioning the lattice £ into A^ sublat- 
tices, such that no neighboring sites are elements of the 
same sublattice. One then chooses any set of A^ mutually 



orthogonal vectors uj^ 
every A^-site simplex. 



e CP^-\ and set z- 
then, each of the uj„ 



= uj On 
vectors will 



FIG. 5: The Q = structure on the Kagome lattice. 



FIG. 6: The V3 X V3 structure on the Kagome lattice. 



occur exactly once, resulting in \Rp\ = 1, which is the 
largest possible value. Thus, the model is unfrustrated, 
in the sense that every simplex F is fully satisfied by the 

assignments, and the energy is the minimum possible 
value: ^ 0. 

For N = 2, there are two equivalent ways of partition- 
ing a bipartite lattice into two sublattices. For N > 2, 
there are in general an infinite number of incquivalcnt 
partitions, all of which have the same ground state en- 
ergy Eq = 0. At finite temperature, though, the free en- 
ergy of these different orderings will in general differ, due 
to the differences in their respective excitation spectra. 
A particular partition may then be selected by entropic 
effects. This phenomenon is known as 'order by disorder' 

To see how entropic effects might select a particular 
partitioning, I derive a nonlinear cr-model, by expanding 
H^^ about a particular zcro-tcmpcrature ordered state. 
Start with 



/■\ /i t \i/2 



(70) 



where ttJ lu 



\Rr 



'(i) 



0. We may now expand 



2 I « cr{j) 

i,3 



(71) 



Thus, the 'low temperature' classical Hamiltonian is 

|2 



LT 



where the sum is over all nearest neighbor pairs on the 
lattice. The full SU(Af) symmetry of the model is of 
course not apparent in eqn. 1721 since it is realized non- 
linearly on the TTj vectors. 

Each TTj vector is subject to a nonholonomic constraint, 

Ti"! ■"'i 5: 1- To solve for the thermodynamics of H^r^, I will 
adopt a simplifying approximation, in which there is just 
one nonholonomic constraint, 7r| tt^ < A/", where hf is 



the number of sites in the lattice. I fix the constraint by 
introducing an auxiliary variable X and demanding 



(73) 



which I enforce with a Lagrange multiplier A. The re- 
sulting model is 



H 



LT 



The local constraints ttJ lo 



(74) 



<t(j) 



are retained. 



It is convenient to rotate to a basis where 
,t 



in which case '^o.(j) tt^ 



r^. ^^.p i.e. the a{i) component 

of the vector TTj. So long as cr(i) ^ a{i) for nearest neigh- 
bors i and j , the local constraints have no effect on the 
Hamiltonian. We arc then left with a Gaussian theory in 
the TTj vectors. Leaving the constraint term aside for the 
moment, we can solve for the spectrum of the first part 
of i^L^p. From this spectrum, we compute the density of 
states per site, g{e). The free energy per site is then 



oo 

^ = -\ + X\x\'' + T jdeg{e) In 



e + \ 
T 



(75) 



Since M is thermodynamically large, we can extremize 
with respect to A, to find the saddle point, yielding 



T di 



gje) 

e + X 



(76) 







Setting A = X 



0, I obtain an equation for T^, 
-I -1 



5(f) 

e 



(77) 



For T < T^, there is Bosc condensation, and |Xp > 0. 
For T > T^, the system is disordered. I stress that the 
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FIG. 7: Free energy for SU(3) simplex solid states on the 
Kagome lattice. Dashed (red): Q = structure. Solid (blue); 
structure. 



disordered phase, described in this way, does not reflect 
the SU(A'^) symmetry which must be present, owing to 
the truncation in eqn. 1721 

A natural setting to investigate the order by disor- 
der mechanism would be the SU(4) SS model on the 
pyrochlore lattice. I defer this analysis, together with 
a companion Monte Carlo simulation, to a later pub- 
lication. Here I will provide a simpler analysis of the 
SU(3) Kagome SS. Since the Mermin- Wagner theorem 
precludes spontaneous breaking of SU(3) in two dimen- 
sions, this analysis at best will reveal which correlations 
should dominate at the local level. The two structures I 
wish to compare are the Q = structure, depicted in fig. 
[5l and the x -s/S structure, depicted in fig. [S] Here, 
the A, B, and C sites correspond to CP^ vectors. 





(78) 



Both the Q = and a/3 x -y/S structures are unfrustrated, 
in the sense that the interactions are fully satisfied on 
every simplex - \Rp\'^ = 1 for all F. Entropic effects, 
however, should favor one of the two configurations. 

The Q = structure may be regarded as a triangular 
Bravais lattice with a three element basis {e.g. a trian- 
gular lattice of up-trianglcs). Each triangular simplex 
contains three tt vectors, each of which has two indepen- 
dent components (neglecting for the moment the global 
constraint). Solving for the spectrum in the absence of 
the constraint, one finds six branches: 



;^^(fe) = 2±2cos(ifc-a 



(79) 



where fc is a vector in the Brillouin zone, and the direct 
lattice vectors are 



a 1,0 



V3^ 



with Og = 02 — fl]^. This results in a free energy per site 
of 



/m = -A(T, + I/<iol„(MI)±l±l£2££,.,81, 





with 



A(r) = V4 + T2 - 2 



(82) 



For the VS X \/3 structure, the underlying lattice is 
again triangular, but now with a nine element basis (see 
fig. [S]). The Hamiltonian iJ^T then purely local, and 
there is no dispersion. The density of states per site is 
found to be 

g{s) = i S{e) + 1 ^(e - 1) + i 6{e - 3) + i 5{e - 4) . (83) 
The free energy per site is 



f{T) = 2-uiT) + lTln{ 



with u{T) = A(T) + 2 satisfying 

1 _ u (u2 _ 3) 



2^6 



T (m2 - 1) (u2 _ 4) 



(84) 



(85) 



Our results are plotted in fig. [T] One finds at all 
temperatures T > that 



Q=oV 



(86) 



(80) 



suggesting that the local correlations should be better 
described by the y/S x \/3 structure. 



VIII. AT THE EDGE 

With periodic boundary conditions applied, the trans- 
lationally invariant AKLT states arc nondegenerate. On 
systems with a boundary, the AKLT models exhibit com- 
pletely free edge states, described by a local spin on 
each edge site which is smaller than the bulk spin S. 
The energy is independent of the edge spin configuration, 
hence there is a ground state entropy {2S^ + 1) InNc, 
where is the number of edge sites. As one moves away 
from the AKLT point in the space of Hamiltonians, the 
degeneracy is lifted and the edge spins interact. 

The existence of weakly interacting S = ^ degrees of 
freedom at the ends of finite S* = 1 Heisenberg chains 
was first discussed by Kennedy p6l |. who found numeri- 
cally an isolated quartet of low energy states for a one- 
parameter family of S* = 1 antiferromagnetic chains. 
These four states are arranged into a singlet and a triplet, 
corresponding to the interaction of two S ~ objects. 
The spin quantum number of the ground state alternates 
with chain length L: singlet for even L and triplet for 
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(a) 



(b) 




FIG. 8: (a) Edge states for the linear chain, (b) (10) edge 
sites for the square lattice, (c) (11) edge sites for the square 
lattice. For the AKLT states, there is an effective free spin of 
length Sc = |M(2 — Zc) on the edge sites (see text). 



odd L. The singlet-triplet splitting was found to decay 
exponentially in L as exp(— L/^), where ^ is the spin- 
spin correlation length. Thus, for long chains, the S = ^ 
objects at the ends are independent. Experimental evi- 
dence for this picture was adduced from ESR studies of 
the compound NENP 

The situation is depicted in fig. [5] for the linear chain 
and for the (10) and (11) edges on the square lattice. 
Recall that each link in the AKLT model supplies one 
Schwinger boson to each of its termini. The spin on any 
site is given by half the total Schwinger boson occupation: 
S = ^{blb^+blb^). Consider first the M = 1 AKLT state 
of eqn. [3] on the linear chain. The bulk sites have total 
boson occupancy n ~ 2, hence S ~ 1, while the end sites 
have n = 1, hence S = ^. If the end sites are also to 
have 5 = 1, they must each receive an extra Schwinger 
boson, of cither spin (| or I). Thus, the end sites are 
described by an effective -S* = ^ degree of freedom. Each 
of these four states is an exact ground state for the AKLT 
Hamiltonian. written as a sum over projection operators 
for total bond spin 5'„ = 2 [H. 

Consider next the square lattice with M ~ 1, for 
which the bulk spin is S* = 2. For a (10) edge, the 
edge sites are threefold coordinated, and each is 'miss- 
ing' one Schwinger boson. The freedom in supplying the 
last Schwinger boson corresponds once again to a S* = i 
object at each edge site. Along the (11) edge, the sites 
are twofold coordinated, and must each accommodate 
two extra bosons, corresponding to S* = 1. The general 
result for the edge spin is clearly 



iA/(z-^J 



where z and are the bulk and edge coordination num- 
bers. I stress that the edge spin configurations arc com- 
pletely degenerate at the AKLT point, since all the inter- 
nal links are satisfied, i.e. annihilated by the local projec- 



tor(s) in the corresponding AKLT Hamiltonian. Moving 
off of the AKLT point, in the direction of the Heisenbcrg 
model, the edge spins will interact. Based in part on 
Kennedy's results, I conclude that the S = ^ chain along 
the (10) edge is antiferromagnetic, while the 5 = 1 chain 
along the (11) edge is ferromagnetic (since consecutive 
edge sites are connected through an odd number of bulk 
sites). 

By deriving and analyzing lattice effects in the spin 
path integral for the Heisenberg model, i.e. tunneling 
processes which have no continuum limit and which do 
not a ppe ar in the effective nonlinear sigma model, Hal- 
dane [28| argued that Heisenberg antiferromagnets with 
25 = mod 4 on the square lattice should have nonde- 
generatc bulk ground states. This result was generalized 
by Read and Sachdev [2§|, who, building on an earlier 
large- Schwinger boson theory [30|, extended Haldane's 
analysis to SU(A^), for {N,N) models on bipartite lat- 
tices. This established a connection to the AKLT states, 
which are nondegenerate in the bulk, and which exist 
only for 5 = 2M on the square lattice. 

For the simplex solids, a corresponding result holds. 
Recall the N = 3, M = 1 model on the Kagome lattice, 
where each site is in the fully symmetric, six-dimensional 
□□ representation, whose wavefunction j^*) is given in 
eqn. [551 Along a (10) edge, as in fig. [21 the edge sites 
each belong to a single simplex. Hence, they are each de- 
ficient by one Schwinger boson. The freedom to supply 
this missing boson on each edge site is equivalent to hav- 
ing a free edge spin in the fundamental representation at 
every site. Thus, as in the SU(2) AKLT case, the bulk 
SU(A^) representation is 'fractionalized' at the edge. The 
□ objects along the edge are of course noninteracting and 
degenerate for the SS projection operator Hamiltonian. 
For general SIJ{N) models which are in some sense close 
to the SS model, these objects will interact. 



IX. CONCLUSIONS 

I have described here a natural generalization of the 
AKLT valence bond solid states for SU(2) quantum anti- 
ferromagnets. The new simplex solid states | ^'(£; M)) 
arc defined by the application of A^-sitc SU(Af) singlet 
creation operators TZ\. to A^-site simplices F oi a par- 
ticular lattice. For each lattice, a hierarchy of SS states 
is defined, parameterized by an integer Af , which is the 
number of singlet operators per simplex. The SS states 
admit a coherent state description in terms of CP^~^ 
variables, and using the coherent states, one finds that 
the equal time correlations in | 'I'(>C ; M)) are equivalent 
to the finite temperature correlations of an associated 



(87) classical CP^ ^ spin Hamiltonian H^^, and on the same 



lattice. The fictitious temperature is T = 1/M, and a 
classical ordering at corresponds to a quantum phase 
transition as a function of the parameter M. This transi- 
tion was investigated using a simple mean field approach. 
I further argued that for N > 2 the ordered structure is 
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selected by an 'order by disorder' mechanism, which in 
the classical model amounts to an entropic favoring of 
one among many degenerate T = structures. I hope to 
report on classical Monte Carlo study of H^^ on Kagome 
{N = 3) and pyrochlore {N = 4) lattices in a future 
publication; there the coherent state formalism derived 
here will be more extensively utilized. Finally, a kind of 
'fractionalization' of the bulk SU(iV) representation at 
the edge was discussed. 
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B. Resolution of identity 



Define the measure 

^ d Re(zj ) d Im(zj ) 



TT 



Next, consider the expression 



y ^ \n)( 

n,---n., 3 



-1 "JV 

oo 



n 



p=0 



where | n) = | rij^, . . . , n^) and 



(92) 



(93) 



r r / ^ \ ^ 

n n .7 — 1 .7 — 1 



XI. APPENDIX I : PROPERTIES OF SU(iV) 
COHERENT STATES 



A. Definition of S\J{N) coherent states 

Consider the fully symmetric representation of SU(A^) 
with p boxes in a single row I call this the p- 
representation, of dimension {^^p~^)- Define the state 

\z;p) = ^^(^z,b\ + ... + z^b^^y\0) , (88) 

where z G CP''^""'^ is a complex unit vector, with z'' z = 1. 
In order to establish some useful properties regarding the 
SU(A^) coherent states, it is convenient to consider the 
unnormalized coherent states 



, ^) = exp z^bl) \ 0) \ z ; p) 



(89) 



Clearly | 2: , ^) is a product of N (unnormalized) coherent 
states for the N Schwinger bosons. One then has 



p=Q 



pi 



Equating the coefficients of {£,£,)^, one obtains the coher- 
ent state overlap. 



(z; p\z' ; p) ^ {z'^z'Y 



(91) 



Here, Xj = \zj\ ■ If we define Xj = (1 — Xj^)yj for j ~ 
1, . . . , — 1, then by integrating over one obtains the 



result 



tN-1 



N-l 



^1 "-N-l 



{N -l + n^ + ... + rij^) ! 



(95) 



Thus, equating the coefficient of in cqn [Ml one 

arrives at the result 



iN-l+p)\ 



p\ 



dn\z; p){z; p \ , (96) 



where the projector onto the p-representation is 

lp= y2 5 \n){n\. (97) 

C. Continuous representation of a state | ?/>) 



Let us define the state 
1 



iibl...,bl)\o) 



^E'^'n&r---ClO), (98) 
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where n = {nj^, . . . , n^}, and where the prime on the 
sum reflects the constraint X^jLi = P- The overlap of 
I ■0) with the coherent state | z ; p) is 



{z\p\'ip) = ^i^nZi' ■■■ ■■ 
n 

= 1p{Zi,...,Zj^) 



■N 



(99) 
(100) 



D. Matrix elements of representation-preserving 
operators 

Next, consider matrix elements of the general 
representation-preserving operator, 



XII. APPENDIX II : THE LOCAL MEAN FIELD 



With the definition of eqn. EH I first compute 



i?^(m) EE Q(iJ A ••• A Q 



(N) 



where 



1 — TO 



S + TO P 



N 



Q/3 — ^fi 



(105) 



(106) 



(101) 



Here the prime on the sum indicates a constraint 
= S.-i = Then 



(102) 



k.l 



N 

n[K ^„ 



i=i 



J J' J J 



where the double prime on the sum indicates constraints 
on each of the sums for m, n, k, and I. This may also 
be computed as an integral over coherent state wavefunc- 
tions: 

/ d^l 4){z) Tj^{z, z) 'ip{z) = ^4>m i ipn 

m.n 

N 

n 



k,l 



m ■ -{-k ■ _n ■ ■ 

X I I Zj ' ' Zj' ' 

N 



k.l 



m-+k- ,n,+i„ 



{<t>\TK\^) 



[N -l+p + K)\ 
Thus, the general matrix element may be written 
(N -l+p + K)\ 



(103) 



d^i(j){z)Tj^{z,z)i/j{z) , 
(104) 



where Tj^{z, z) is obtained from eqn. llOll bv substitution 
6,- 2„ and bl — * 



is the projector onto subspace spanned by w'^. I now 
systematically expand in powers of the projectors and 
contract over all free indices. The result is 



i?^(TO) = iV! ^ — 



N 



The local mean field on a j ~ 1 site is given by the 
expression in eqn. 1611 Expanding the numerator. 



"jv/3]v "iv-i/3jv-i 

(107) 

in powers of the projectors, the term of order j is 



{N-j-iy. 



1 — TO 

N 



N-j-l 



/c^<...<A:j.< 



(108) 



TV 



Writing 



/3,/32.../3j, 



= > sgn((T) 5 



(109) 

we see that once this is inserted into eqn. 11081 the only 
surviving permutations are the identity, and the {N — 1) 
two-cycles which include index N . All other permuta- 
tions result in contractions of indices among orthogonal 
projectors, and hence yield zero. Furthermore using com- 
pleteness, we have 



N-l 

1=1 



(110) 
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Thus, From this expression, I obtain the results of cqns. 

[Ml and EH 



ga^a2...Q„ g/3i/32.../3„ p^i pfc2 ...p""-! 

fci<...<fcj<W 



k^<...<k^<N '=1 
^^-'V ■ (111) 
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